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unifying approach to the theory of random matrices. Introducing an appro- 
priate time t-dependence in the probability distribution of the matrix ensem- 
ble, leads to vertex operator expressions for the n-point correlation functions 
(probabilities of n eigenvalues in infinitesimal intervals) and the correspond- 
ing Fredholm determinants (probabilities of no eigenvalue in a Borel subset 
E); the latter probability is a ratio of r-functions for the KP-equation, whose 
numerator satisfy partial differential equations, which decouple into the sum 
of two parts: a Virasoro-like part depending on time only and a Vect(S' 1 )-part 
depending on the boundary points A; L of E. Upon setting t — 0, and using 
the KP-hierarchy to eliminate t-derivatives, these PDE's lead to a hierarchy 
of non-linear PDE's, purely in terms of the Ai. These PDE's are nothing 
else but the KP hierarchy for which the t-partials, viewed as commuting op- 
erators, are replaced by non-commuting operators in the endpoints Ai of the 
E under consideration. When the boundary of E consists of one point and 
for the known kernels, one recovers the Painleve equations, found in [14,24]; 
from the present work, it also appears that some of the Painleve equations 
can be viewed as the KP equation in non- commutative operators. Some of 
these results have been announced by us in [5]. 

Random matrices provide a model for excitation spectra of heavy nuclei 
at high excitations (Wigner [27], Dyson [12,13] and Mehta [18,19]), based on 
the nuclear experimental data by Porter and Rosenzweig [21]; they observed 
that the occurrence of two levels, close to each other, is a rare event (level 
repulsion), showing that the spacing is not Poissonian, as one might expect 
from a naive point of view; this lead Wigner to his so-called surmise. 

In their pioneering work, Jimbo, Miwa, Mori and Sato[14] have shown 
some (p, g)-variables derived from the sine kernel satisfy a certain Neumann- 
like completely integrable finite-dimensional Hamiltonian system. They also 
showed that the distribution of the level-spacings for the sine kernel satis- 
fies a Painleve equation. Tracy and Widom [24,25] have successfully used 
functional-theoretical tools to compute the level spacing distributions for a 
more general class of kernels, always yielding Neumann- type systems. Deift, 
Its and Zhou [10] have used the Riemann-Hilbert approach to find the pre- 
cise asymptotics for the distributions above. Random matrices have come up 
in the context of statistical mechanics and quantum gravity; see [7, 15, 17]. 
Random matrix ideas play an increasingly prominent role in mathematics: 
not only have they come up in the spacings of the zeroes of the Riemann zeta 
function, but their relevance has been observed in the chaotic Sinai billiard 
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and, more generally, in chaotic geodesic flows; Sarnak p4[| conjectures that 
chaos leads to the "spectral rigidity" , typical of the spectral distributions of 
random matrices, whereas the spectrum of an integrable system is random 
(Poisson)! 

As is well known in the theory of random matrices, probability distribu- 
tions of finite matrix ensembles are usually given by explicit matrix integrals, 
whereas the probabilility measure for infinite matrix ensembles is given by 
certain kernels. We found out that the first situation, actually leading to 
Christoffel-Darboux kernels, typically ties up with the Toda lattice and the 
second, corresponding to kernels represented by an integral, ties up with the 
KdV equation. This paper deals with the second case (KdV) only. There 
is a striking analogy between the two cases, showing our point of view is 
very robust. A similar theory can be worked out for finite matrix ensembles 
(discrete kernels and 1-Toda lattice) and for two coupled random matrices 
(2- Toda lattice). Theorems 0.1 and 0.2, stated in this section, are the main 
results of this paper. 

In this work, we shall be dealing with the KP hierarchy, briefly ex- 
plained in section 2; it is a hierarchy of isospectral deformations of a pseudo- 
differential operator L = D + Yli>i a i{ x ,t)D~ l , with[| D := d/dx, 

BT 

— = [(L») + ,L],for teC 00 . (1) 

We also consider the p-reduced KP hierarchy, i.e. the reduction to L's such 
that L p is a differential operator for some fixed p > 2. Sato tells us that 
the solution L to equations (1) can ultimately be expressed in terms of a 
r-function. The wave and adjoint wave functions, expressed in terms of the 
r-function[] 

¥(*,*,*) = e- + ^^ i T(t and 9*{x,t,z) = e~"~^ ^ + }*~ ' D ; 

v 1 r(t) v ' ' 1 r(t) 

(2) 

satisfy 



2 A + denotes the differential part of the pseudo-differential operator A 
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In this paper, as in the general theory of integrable systems, vertex op- 
erators play a prominent role : they are Darboux transforms involving all 
times. In particular, for the KP equation, the vertex operator 

X(t, y, z) := ^ e ^-ynu e ^T(y-'-^A (4) 

has the striking feature that X(t,y,z)r and r + X(t,y,z)r are both r- 
functions. Given distinct roots of unity u,u' G ( p := {u \ u p = 1}, the 
vertex operator X(t,uz,u'z) maps the space of p-reduced r-functions into 
itself; for vertex operators and their interaction with the Virasoro algebra, 
see sections 3 and 4. 

Sections 3,4, 5 and 6 will lead to the first main theorem, the first part 
of which is a continuous version of a "soliton" formula. These sections and 
Theorem 0.1 are due to Adler-Shiota-van Moerbeke ||. For the sake of 
completeness, we sketch the proof in this paper. 

Theorem 0.1. Define the (x,t)-dependent kernel k Xjt {y,z) and k^ t (y,z) : = 
k x j(y, z)I E (z) with x G E, t G C°°, y, z G C, and E C R + a Borel subset: 

/X 
dx^2a w ty*(x,t,ujy) ^ b^{x, t, u'z), (5) 

where ^(x,t,z) and ^*(x,t,z) are the wave and adjoint wave function for 
the p-reduction of KP and where the coefficients a u , G C are subjected to 
J2u>e( P ^Or = °- Then the following holds: 

(i) the kernel k(y,z), its determinant and its Fredholm determinant are all 
three expressible in terms of the vertex operator 

Y(t,y,z):= 22 aubu>X(t,uy,uj'z) (6) 
acting on the underlying T-function^ : 

T 

4 Y(i, z, z) in formula (7) below must be understood as lim y ^ z Y(t, y, z). 
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det(k Xti (y i} Zj))i<ij< n = -Y[t, y 1} zi) . . . Y(t, y n , z n )r 

det(J - nk* t ) = ^- e -^z dz Y ^ z \ (7) 

("continous" soliton formula) 

(ii) Let the kernel k x ,t(y, z) in (5) be such that the underlying r-function of 
\l/ and ty* satisfies a Virasoro constraint 

WtyT = CkpT for a fixed k > — 1. 

Then for the disjoint union E = U[=i[ a 2i-i) a 2i] C M + , the Fredholm deter- 
minant det(I — fik^ t ) satishes the following Virasoro constraint for that same 
k > -1: 



2r f) 

E a " P+1 ^ + - c ^ ) (rdet(J - /< t )) = 0; (8) 

i=i * / 

note the boundary a-part and the time t-part decouple. 

In the next theorem we apply equation (8) to compute the partial dif- 
ferential equations for the distribution of the spectrum for matrix ensembles 
whose probability is given by a kernel. To state the problem, consider a 
first-order differential operator A in z of the form 



a = a z = \z-^ (|+n*))+l>_ 2i 



with 



(9) 



a 8 

V(z) = -z + ^z 3 ^0, m = deg V = or 2, (10) 
2 6 



5 Define W 7 ? ] = <L 



d/dtn if n > 



if n = i+j=n 



4 X) :=^ X) = ifn<0 , j( 2 > := + („ + l)W^ = £ ^ 
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and the differential part of its "Fourier" transform 

d 



A = A X = (^(x + V'(D))D- m+1 + J2^2iD~ 2i \ with D — 

(11) 

r 

Given a disjoint union E = \^J[A 2 i-i, A 2 i] C R + , define differential operators 

i=l 

A n , which we declare to be of homogeneous "weight" n, as follows 

2r 

E n + l-m (j 
A l 2 — , n= 1,3,5,.... 

We now state the second main theorem, established in section 8: 

Theorem 0.2. Let ^(x, z), x G R, z G C be a solution of the linear partial 
differential equation 

A g *{x,z) = A x *(x 1 z) 1 (12) 

with holomorphic (in z~ l ) initial condition at x = 0, subjected to the follow- 
ing differential equation for some a,b,c G C, 

(aA 2 z + bA z + c)*(0, 2) = z 2 ^(0, z), with (0, 2) = 1 + O^" 1 ). (13) 

Tien 

• ^(x, z) is a solution of a second order problem for some potential q(x) 

(D 2 + q(x))^(x,z) = z 2 ^(x,z). (14) 

• Given the kernel 

K x (y,z) :=I E {z) J 2j/ i/4^i/4 

with 

$( X)M ) ;= b u e uV ^{x,UV,), 
w=±l 

the Fredholm determinant f{A\, y4 2r ) := det(J— XK X ) satisfies a hierarchy 
of bilinear partial differential equations f] in the for odd n > 3: 

i+j=n+l 

6 the Pi are the elementary Schur polynomials e^^° tiZ ' — p n (t)z n , and j?i(±.4) := 
p 4 (±A,0,±^ 3 ,0,...) 
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+ (terms of lower weight i for 1 < % < n) = 0, (16) 
where x appears in the coefficients of the lower weight terms only. 



1 Random matrices reviewed 

Define on the ensemble TCn = {N x N Hermitian matrices} the probability [] 

P{M e dM) = ce- TvV(M) dM, 

where c is a normalization constant. Then for z\, . . . , z^ e R we have Q 

P(one eigenvalue in each [zi, Zi + dZi), i — 1, . . . , N) 
= cvol{U{N))e-^ v(z ^A 2 {z)dz 1 ■■■dz N 

and for < k < N 

P(one eigenvalue in each [zi, z% + dzi], i = 1, . . . , k) 

= c' I j ■ ■ ■ j ' v{Zl) A 2 (z)dz k+1 - ■ ■ dz N \ dzx---dz k 

\ K N - k / 

= c" det(K N (z h z j )) 1 ^ ij ^ k dzi- ■ -dz k , (17) 
and if J C R, then 

(— 1)* d k 

P(exactly k eigenvalues in J) = — — — det (/ — XK^j | A _ 1 , 

where 

K J N (z,z')=K N (z,z')Ij(z>), 

Ij the indicator function of J, and is the Schwartz kernel of the orthog- 
onal projector C[z] — > C + Cz + ■ ■ ■ + <Cz N ~ l with respect to the measure 
e~ TlV ^dz, namely 

JV-l 



K N (z,z') = ^(p k (z)(p k {z') 



k=0 



i<i<j<N dRe(Mij)dlm(Mij) 
8 A(z) = Yli<i<j<N( z i ~ z j)i ^ ne Vandermonde determinant 
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in terms of orthonormal functions ipk{z) = e TrV ^ z '' 2 Pk{z) with respect 
to dz or orthogonal polynomials Pk{z) = (l/y/hk)z k + ••• with respect to 
e~ TTV ^dz. 

In the last equality (=) in (15) one has used the simple property of a 
Vandermonde that 

A(x) 2 = det(? fe _i(xi))i< i)fe < JV 

N 



det ^2lk-i(xi)q k -i(xj) 



l<i,j<N 



in terms of any monic polynomials qk(x) of degree k. Then in (=) the or- 
thonormality of (p k is used to obtain 

K N (x, y)K N (y, z)dy = K N (x, z), 

which just means that gives a projector. When Viz) is quadratic, we 
have for large N, 

P(an eigenvalue G [z, z + dz]) = K N (z, z)dz 

\{2N - z 2 fl 2 dz if \z\ < (2iV) 1/2 
if \z\ > {2N) 1 ' 2 

is given by the circular distribution (Wigner's semi-circle law), we have that 
for z ~ the average spacing between the eigenvalues near the origin is 
~ (Kn^Q))- 1 = n/V2N and near the edge (z ~ V2N) is 1/(2 1 / 2 N 1 / 6 ), 
leading to 

1 / z z' \ , lsm7c{z-z') 

lim ; r-fVjv ; r, : r = K (Z, Z = 

^^(0,0) N \K N {0,0)' K N (0,0)J v ' ; 7t z-z' 

(bulk scaling limit) 

1 / z z' \ f°° 

& Kjp)*'" + ^ + J = I A( X+ z)M X+ z')d X 

(edge scaling limit) 

in terms of the classical Airy function. In a similar context one also finds 
the Bessel kernel; for background on such matters, consult Mehta's excellent 
book 0. 
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2 KP revisited 

Consider the infinite-dimensional span 

W = spanJ^O?) = z k (l + Oiz- 1 )), k = 0,1, 2, ...} 

of meromorphic functions ipk{ z ) i n z ' 1 with poles of order k at z — oo. With 
Sato, define the r-function 

r(x,t) = r(t) := det Proj (V E ^ 2 ' -> H + ) 

where t — t + (x, 0, 0, ...) = (x + £i, t 2 , ...)• The r-function is shown to satisfy 
the following bilinear relation, 

j e £rft-«i>V(t _ [z- 1 ])-^' + [z" 1 ])^ = 0, (18) 

where the integral is taken over a small circle around z = oo; this relation 
actually characterizes the r-function. 

The tau-function r(t) above leads a wave operator 

n=0 

and also to a wave and dual wave function 

*(x,t,z) := (Se**)^"^ = e^+Er^ lff ~ [ f ' ]) , (20) 

r(t) 



r(t) 



(21) 



From (18) it follows that, along a small contour z = oo, 

/ ty{x,t,z)ty*(x',t',z)dz = 0. (22) 

»/ 2 = 0O 



3 i?+ = span c {l,z,z 2 , ...} 
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Conversely, if functions ^ and with the appropriate asymptotic behav- 
ior, satisfy the bilinear identity (22), then \1/ is a wave function and 1 i r * an 
adjoint wave function of the KP hierarchy; i.e., there exists a monic first 
order pseudodifferential operator L = D + • • • = S D S* -1 in x, depending on 
t, satisfying (1) and (3). For more details on these matters, see the work of 
Date, Jimbo, Kashiwara, Miwa and for vertex operators and symmetries, 



Since x and t% always appear together as x + ti in \l/ and in every 
formula involving ^ or and r, we always have t instead of t in the argu- 
ment of r. To simplify the notation in such a formula, we will often omit 
underlining and simply write t instead of t in the argument of r, if there is 
no fear of confusion. 

We also introduce the Orlov-Schulman pseudo-differential operator M, 



satisfying [L, M] = 1, and 

— = Mm. 

oz 

The p-reduction of KP (p > 2) is one where LP is a differential operator; 
it is equivalent to the property that 



see [1,0. 



oo 




(23) 



i 




satisfies z p W C W. 



For future use, we also introduce the 5-function, 



oo 



1 



t 



-1 




n=— oo 



with the customary property 



f(\,l*)6(-) = f(\,\)6(-). 




(24) 
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Note the function 



has the usual ^-function property 

^- J f(X^)5(X,^ = f(X,X) 
and is a function of A — /i only, since 

The action of D~ x on "oscillating functions" like ^* is usually defined 
formally by D~ x e xz = z~ 1 e xz and the Leibniz rule. E.g., if f(x,t,z) behaves 
nicely near z = oo, then 



l—k „xz 



D- 1 (f(x,t,z)e xz ) = jr( k) Dkf ' D 

oo 

= ^2(-l) k D k f ■ z~ l - k e xz . 

k=0 

When this interpretation is not sufficient, e.g., if z needs to be finite rather 
than defined around oo, it is defined so that it gives back this formal inter- 
pretation asymptotically as z — > oo, e.g., 



D- 1 (f(x,t,z)e xz ) 



/ f(x,t,z)e xz dx ifRez>0, 

J — oo 

POO 

- f(x,t,z)e xz dx ifRez<0, 

J X 

/X 
^(x,t,/j,)^*(x,t,X)dx ifRe/i>ReA, 
-oo 
poo 

- V(x,t,[j,)ty*(x,t,\)dx ifRe/i<ReA, 

J X 

if the right hand side makes sense. Hence we shall also denote D' 1 by J dx 

when it is acting on an oscillating function. The lower limit of this integral 
sign, while not written out, is important and must depend, e.g., on the sign 
of Re z or Re \i — Re A in these examples. 



D-\V(x,t,n)V*(x,t,X)) = 
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3 Vertex operators for the KP equation 

Define the vertex operator as in (4): and more generally 
X(t, Ai, . . . , A n ; H\i ■ ■ ■ i AO 



-.YlXiX,,^): (26) 



n n 

with the customary expansion of X(t, A,/x) in terms of iy-generators: 

X(t, A, „) = -L- Y: £ \-"W?\ (27) 

Observe that and coincide with the generators in footnote 3. 

Vertex operators X(t,X,ji) are known to satisfy the following commutation 
relations, due to [[10] : 



[X{\, n),X{u, v)\ = -5{u, /i)X(A, v) + 8{X, v)X{u, fj). (28) 

We also define an operator N, which has a formal expansion similar to the 
vertex operator expansion (E7), but with replaced by z lJrk ~ x [d/dz) k 1 : 

n{X,fi,z) := 5{\,z)e^- x) & 



^ ^^(/i-A)^ 1 ( d 



- — OO 



A* ^ 



^ k=0 l=-oo v 

1 sr^ iPjz ^ k ST \-e-k h j+k-i ( 9 V 

z- — n t— — , \ ' 



fc-1 

(29) 



u — A 

r k=0 

The next Lemma shows the precise relationship between the vertex op- 
erator (|2l]D and the composition of vertex operators (4): 

10 where : : denotes normal ordering, i.e., always pull differentiation to the right. 
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Lemma 3.1. If, possibly after some relabeling, \yi\, \zi\ < 1 2/2 1 , \z 2 \ < 
then the following relation holds: 

n(n — l) 

X(yi;z 1 )...X(y n ;z n ) = (-1) 2 A(z)A(y)X(y 1 , ... 

1 Vn 1 *1 j • ■ ■ j Z n ). (30) 

Proof : Whenever the inequalities above are satisfied, the vertex operators 
automatically commute. Therefore it suffices to prove ([30]) for \y\\, \z\\ > 
\y 2 \, \z 2 \ > setting 

1 

c(y,z) :-- 



(zi-yi)lH,e=2( z k-ytY 
we proceed by induction: 

X(y u z 1 )X(y 2 , ...,y n ; z 2 , z n ) 



c(y,z)e 



n 



-!) n i n^ i_ ^)^ i_2:fc ^^ i '-' ?/n '^ 1 '-'^ 



fe=2 

where we have used 



e ^1 % = 1 — a \a\ < 1 and | — , — , — , — < 1; (31) 

Vi zi Vi z\ 

this ends the proof of Lemma 3.1. 
Corollary 3.1.1. Away from y = z, we have 

X(y,z)X(y,z) = (32) 

and so 



?aX(y,z) 



l + aX(y,z). (33) 
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Proof : Indeed from (|30|), we have, assuming \yi\, \zi\ < \y2\, \z 2 \ 

X(y 1 ;z 1 )X(y 2 ;z 2 ) = -(yi ~ z 2 )X(y 1 ,y 2 ; Zi, z 2 ) 

and thus, in the limit, when y\ — > y 2 , z\ — > z 2 , we find, away from 
yi = zi, the vanishing (|32|) of the square of the vertex operator. Thus Taylor 
expanding e aX{ ^)^ leads to (|33|) . 

Lemma 3.2. jj^ The following commutation relation holds: 



v t + kldY_(_d\ k ou l+ky X{lir] 



dv J \ du J / ' k + 1 



Wi k+1 \x(u,v) 



(34) 

Proof : First observe, using the property (|24] ) of the ^-function, that 



n{\yL,v)X{u,v) = 8(X,v)e^- x) &X(u,v) 
= 5(X, v )X(u, v + n — A) 
= 5(\, v )X(u, n) 

and, using 5(y, z) is a function of y — z (see (|23|)), that 



n(X,u,,u) T X(u,v) = e (A - M) ^<5(A,w)X(w,t>) 

= <5(A, m + A — \i)X{u + X — fi,v) 
= 5(u, n)X{u + A — fi,v) 
= 5(u, /x)X(A, v ). 

From the commutation relation (p8| ) and the the two relations above, it 
follows that 

[X(X,fi),X(u,v)] = -5(u,fi)X(X,v) + 5(X,v)X(u,fi) 

= (-n(X, fi,u) T + n(X,fM, v)) X(u,v). 

Finally, expanding both sides of 

[X(X, n),X(u, v)\ = (n(A, fj,, v) — n(A, //, u) T )X(u, v) 
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in A and /i according to fl2"T|) and (29), yields 



E 

k=l 



fl — X 

(m-a)*- 1 " 



E 

.fe=0 



(a* -A)' 



A:! 



k\ 



X~*~ k k 



\- £ ~ k wi k \x( 

e=-oo 



9m' 



which, upon comparing coefficients, ends the proof of Lemma 3.2. 

Corollary 3.2.1. For all k £ Z, we have 

(v e - v) = [wf \ X(u, v) 

and 



Wi 2 \x(u,v) 



(35) 



(36) 



which further translates into the following statement for "(1/2, 1/2) -differentials" , 
(v e - u e )X(u 1 v)Vdudv = wj; 1 \x(u } v)Vdudv 



(u e+1 — + v e+1 —)X(u, v)Vd^h = Wi 2) + (I + l)wl l \ X(u, v)Vd^h 
ou ov L 



Proof : The first two formulas follow at once from setting k = and 1 in (|34|) . 
The third relation follows immediately from (j35|). The last one follows from 
combining the two computations below: on the one hand, both (|35|) and (j3G|) 
lead to 



jP,X(u,v) 



1 



Wjf ] + (k + l)W^>,X(u,v 



(i) 



k+l9 ' 9 - k+1 ' k+1 (v k -u k ))X( Uj v) 



dv ^ du U 2 
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on the other hand, we have 



fe+i 



d_ 

du 



){f{u)Vdki) 



—f(u + eu k+l Wd(u + eu k +^) | £=0 
^ (/(«) + eu k+1 f'(u)) (l + £ -{k + l)u fc ) 




e=0 



and similarly for u replaced by v. Setting / = X(u,v) ends the proof of 
Corollary 3.2.1. 

4 Virasoro equations satisfied by vertex op- 
erators and r-functions 

For distinct roots G ( p , the vertex operator X(t,wz,u'z) transforms 

the space of r-functions for the p-reduced KP into itself and satisfies some 
infinite set of differential equations, involving the Virasoro algebra; for the 
definition of Wf>, see footnote 3. This section is due to Adler-Shiota-van 
Moerbeke §. 

Theorem 4.1. For distinct £ ( p and for i e Z, i > —p such that p\l, 
the vertex operator Y(z) := X(t,uz,u'z) satisfies the equations: 




e -\f*dzY(z) 



dzY(z) 



0. 



0. 



(37) 



(39) 



(38) 
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Proof : At first we compute, using Corollary 3.2.1, uj p = u' p = 1 and p\£: 



^(z^X^u'z)) 
z e+1 -^ + (£+l)z e ^j X(uz,u'z) 

z e+1 co^- + z i+ v£ + (£ + l)z e ) X ( ,,. r ) 
ou ov 

u e+1 ^- + v e+1 ^- + (£ + iy ) X(u,r)\ 



V — UJ f Z 



UJZ 

v—uj'z 



wi 2) ,X(u,v) 
Wf\X{uz,dz) 



U — UJZ 
V = UJ f z 



1 

2 L 



U — LUZ 
V — ijj'z 



establishing fl3"7j). Upon using the above relation, one computes 



db da + 2 



dzX(uJz, lu'z) 



wf\x{uz,dz) 



= -b i+1 X{ujb,u'b) + a i+1 X{uja,u'a)+ I dz 

J a 

= -b e+1 X(ub, u'b) + a i+1 X{ua, uo'a) + z e+1 X{uz, J z) \ b a = 0, 
establishing (|3~8l). Note 



V e :=-b^^--a i+1 
ob 



1 + ^ 



is a derivation. Therefore upon setting U := J b Y(z)dz, expanding the 
ponential and using V(U = 0, we have 



V £ (e 



III 



J2 J2 u n - k {v t u)u k - 1 = 0, 



k=l 



establishing (39) and Theorem 4.1. 
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Corollary 4.1.1. The vertex operator 

Y(z) := e v ^X(t,uz,u'z) 
with distinct u, u) 1 £ Q p satisfies the equations: 



d_ 

dz 



z e+1 f(z)Y(z) = Wf>,Y(z) 



(2) 



W, 



(2) 



(2) 



(40) 



dzY(z) = 0. (41) 



where 



W 



(2) 



fe>0 



J V 0^]^+— V Wj+Li' for £ £ Z, £ > -p suci tiat p\£ 

fe>0 



V'=: 



/(*) E t >o a p^ pi ' 



Proof: 



^-z i+1 f(z)e v ^X(uJz,uj , z) 
^(z)/^ / +1 X(^,^)e y ^ +e v ( 2 ^(/(z)^ +1 X(^z,^)) 



b k z k+E+1 X(uz,u'z) J + 



a k z k+e+1 X(uz,u'z) 



UJ — UJ 



' k+e+i ' 



X(uz, u'z)e 



kepZ 

(by ([35|) and 



uj k+e+1 = UJ 



UJ 



tk+e+i _ , ,/ 



u/) 



+- 



J2 a M% , X(^,^z)e y(2) 



(by (BS 
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concluding the proof of relation ( EDI ) while the two relations following (ffl|) 
and ([42|) are proved exactly in the way the relations (^8|) and (^) are derived 
from (0). 

The next theorem shows that if a function / satisfies a Virasoro con- 
straint, then the Virasoro vector fields acting on the function 

„— A f b dz X(ujz,lu' z) r 

move the end points a and b according to the simple vector field 

a = a kp+1 and b = b kp+1 . (43) 
Theorem 4.2. For an arbitrary function f(t) satisfying for some k > — 1 

Wgf = c kp f, 

we have for that same k 



db 



Proof : Setting, as before, U = J dzX(tuz, uj'z), it follows from Theorem 4.1 
that 







db 



da 2 



wi 2) 

yv kp ) 



e- XU f(t) 



f_ b kp+l^_ _ Jp+ll. , l w (2) _ 1 A e ~XU f(t) 



using in the last equality the hypothesis, thus establishing Theorem 4.2. 



5 Fay identities and symmetries 

The bilinear identity ( fl8|) generates many identities for the r-function. In 
the next proposition we describe some of them, which will be useful in the 
next section: 
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Proposition 5.1. (Fay identities). The T-function satisfies the following 
bilinear identity 

( s ° - s i)( s 2 - s 3 )r(t + [s ] + [si])r(t + [s 2 ] + [s 3 ]) = 0. (44) 

cyclic permutations 
of indices {1,2,3} 

and a more general identity 

'<t + [y- '] - [zj l \) 1 



det 



( T(t) Hi - Zj J l<i,j<n 

n(n-D A(y)A{z) r ( t + Eik rl ] -EiK 71 



(45) 



Proof : Upon shifting t appropriately and using the residue theorem, the bi- 
linear identity (18) leads to the quadratic Fay identity It is also the simplest 
case of the Pliicker relations for the corresponding infinite-dimensional plane 
m, |2§ 0). 



sec 



The proof of the second relation ( P5| ) proceeds by induction and is due to 
; at first, it is based on the bilinear relation (HI) for 



t^t+fo 1 } and t'^t + J^y- 1 ]- J>7 X ]- 

k=2 l=\ 

Using (PH), one computes 

// n n \ 

V 2 £=1 / 



// n n \ 

\ fc=2 <=1 / 

ie- (i - ; 



nu (i - 



-dz 
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-r(t)r t + E^-E^ 



a vi u: = 2 (i 



ML 

Ha 



k=l 



7=1 



k=2 



- if) n « 



(46) 



Develop now the determinant on the left hand side of (f45|) according to 



the first row and assume the validity of (J5) for n i— > n — 1; we find 

y, yt- At + fa 1 ] - [zj 1 ]) 1 ^ / T (t + [yr 1 ] - K 1 ]) 1 
r(t) Vx-Zj \ 



T(t) 



Vi - z k 



2<i<n 
1 < k < n 



E, 1 y 7 -i r ( t + [Vi *] ~ l z j *]) 1 / (n-iKn-2) A(y 2 , -,1/n) A (2:1, 4, ...,2: n 



3=1 



Vi-z 



n^ 1 (yk 



rit + E^yA-Ek^k 1 '. 



it) 



A(y)A(*) r^ + Eik^-Eit^ 1 ] 



using the identity (|46|) ; note for n = 1, relation ( |45f) reduces to ( B) , ending 
the proof of Proposition 5.1. 

If r is a r-function, then X(y, z)t is also a r-function, but also 



D aX(y,z) 



t = t + aX(y, z)t, 



as follows from fl33|). But more is true, and it is an important feature of 
the vertex operator: X(y, z) is a generating function of the infinitesimal 
symmetries of the manifold of r-functions ]T0| . 

A generating function for the symmetry vector fields on the manifold of 
wave functions ^(x, t, z) is given by 

Y N V = -NJil , with N(z,y) := e {z ~ y)M 5 (y, L) . 
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where N is the pseudo-differential operator corresponding to n, defined in 
(29), by means of the relation n\l/ = N^f; i.e., z and d/dz in n get replaced, 
in reverse order by L and M defined in (f^J). Note the operator identity 

N(z, y)_ = -V(x, t, z) o D- 1 o i&*(x, t, y). (47) 

Upon using the shift operator^] r] = T](z) = ^^(l/i),2~ l <9/<9tj, the relation 
between the infinitesimal symmetries represented by X on the space of r- 
functions and the ones represented by Yn on the space of wave functions is 
given by (see || and, for generalizations, [pH]). 

Lemma 5.2. (ASV identity). 

where * = z), r = r(t), iV = JV(tt, y) and X = X(y, u). 



6 Determinant of kernels and vertex opera- 
tors 

The main contention of this section is Theorem 6.1; at first, it claims that 
the kernel 

D- l (^*(x,t,y)^(x,t,z)) 

with D^ 1 interpreted as in section 2, is the vertex operator applied to the 
r-function divided by r. Secondly, the determinant of a matrix involving this 
kernel equals the vertex operator iterated several times as in fl30f ) . Theorems 
6.1 and 6.2 are due to ||. 

Theorem 6.1. (KP) 
(i) 

D- l (**(x,t,y)*(x,t,z)) = -±-X(t,y,z)r(t). (49) 
T { t ) 

11 e-V{t) = f{t-[z]) 
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(ii) Assuming \yi\,\zi\ < 1 3/2 [ , | ^2 1 < ••• < \y n \, \ z n\, we have the following 
vertex operators relation 

d^(D- 1 ^{x,t,y i Mx,t,z j )) = det(-^X(£, M )T(0) (50) 

" \ T K l ) J l<i,j<n 

= - X(t,y 1 ,z 1 )...X(t,y n ,z n )r 

(51) 

Proof : The identity in (i) established in |I| can be obtained at once by taking 
a limit of identity (0), when u — > z. Indeed, from (f47|) and the prior identity, 
we have for the left hand side of (BSD, 



Y N V(x,t z) = _ D -iy*( t ^ t z) 



U^rZ 



and for the right hand side Q, one finds, upon using the explicit formula (4) 
for the vertex operator X(t,y,u) and (|3l|) 

limfc"" i\ X ( t >y> u ) T ( t ) = i im ( z -uX(t,y,u)e- r >r X(t,y,u)r(t) 



r(t) u^zyz — y e V r(t) 

X{t,y,z)r{t) 
r(t) ' 

thus leading to the desired identity. 

Another proof is based on Fay's trisecant identity (f44|); indeed differen- 
tiating ( ^4|) with regard to so, setting so = S3 = 0, dividing by sis 2 and 
shifting t —*t — [s 2 ], yields 

{r(t),T(t+[ Sl }-[s 2 })} 

+(s^ - s^ 1 ) (r(t + [ ai ] - [s 2 ])r(t) - r(t + [ Sl })r(t - [s 2 ])) = 0. 

Setting y = sj" 1 , z = s^ 1 and multiplying with r(t)~ 2 exp ^ ti{ z% — y l ) lead 
to 

T(t + [y-i]) c _ y ^ t ^ T(t-[z-i]) cyruzJ = ( ^^ rjt + [y- 1 ] - [■ 

r(t) r(t) z-ydx \ r(t) 

12 use rt/§J(t^j = *E%X(t, y, u) 
13 {/, 5} = Wronskian[ 5 , /] = f'g - fg'. 
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which implies (i). 

To prove (]5lD , multiply the higher Fay identity (45) with appropriate 
exponentials, which leads to the following vertex operators relation: 

det h^^^w) = 2 A 0) A u/) 



l<i,j<n 

1 



- X(t,yi;zx)...X(t,y n ;z n )T, 



T 



using in the last equality the identity (PU|), which is valid under the assump- 
tion | yi |, \zi | < 1 2/2 1, 1 32 1 < upon using (|49|) in the left hand side of the 
above identity, one finds flSTD ,ending the proof of Theorem 6.1. 

For future applications, we consider more general kernels of the form 
D- 1 aj$*(x,t,u)y) ^ b^(x,t,cu f z), 

involving coefficients a u and b w G C, uj G ( p , satisfying 

E a wbuj _ 
to ~ 

This kernel can again be expressed in terms of a basic vertex operator 
Y(y,z):= ^2 aJ)u>X{t,uy,u'z) 

leading to a theorem, completely analogous to Theorem 6.1. 

Theorem 6.2. (p-reduced KP) 

(i) 

D- 1 a u **(x, t, toy) £ b w*(x, *, J z) = ±-Y(y, z)r. (52) 
(ii) If 1 2/1 1 , |*i | < \y 2 \, \z 2 \ < then 
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det I D^C^a^it.uy^Y^ b^(t,u' Zj )) J = det Y( Vi , Zj )r 



l<i,j<n 

Zi) 



r ( i Y(V> 

V ordered 



If 

^ ^ = 0, (54) 



the right hand side of ( |53|) has no singularities in the positive quadrant {yi > 



and Zj > with i,j = 1, n} and lim^^ Y(y, z) exists. 

Proof of Theorem 6.2 : The proof of (i) follows at once from Theorem 6.1. 
To prove (ii), we set 

k(y,z) := X(t '^ Mt) = D-^*(t,yMt,z); 

we have 14 

det lD-\Y^^*^yi)){^2 b^(t,Lu' Zj )) 
= det( ^2 oJ>u'k(uyiitJz j ) 

n 

= E(- i ) £(<r) n( e ^k^jz^) 

(T(zTr n i=l uj,u}'£^ p 

(-l) e(<T) II a^b^k^yU^h^) 

—t—ff, . Ki<n 



l<i,j<n 



l<i ,j<n 



l UJkU = (Wfelj/l, OJknUn) 
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upon interchanging product J^J and sum 

i=l u) ,u>' &C, V 

— »' — >' 

upon replacing u e (Cp) n bycr co> e (Cp) n ; 

— >/ — >/ 
note (<TU))i= {uj ) a . 



-j, 



u,u e(Cp) r 



l<i,j<n 



ri n 

e n a ^ n V'),- det ( A; (^)*2/*' 



l<i j'<n 



j n n n 

7 e nv)iiiV)i(n ^.(^ii^ 

W,W e(Cp)" ordered 

for 1 , I < 1 2/2 1 , W\ < using \(u)i\ = 1 
1 n 

~( TT ( ^ a^X^,^))^ 



r 

ordered 

1 



Y(y h z 1 )Y(y 2 , z 2 )...Y{y n , z n )r 

T 



ending the proof of fl53|). 



We now prove that the expression above 

rdet IT 1 ^a w tf b u ^(t,u' 



l<ij'<n 



e ri" ' n x ((£i)iyi, (w 2 )»zi) j r 

iZi,<3 2 e(C P ) n i=l J'=l \ «=i / 

\ ordered / 

n(n-i) TT TT 7 A(w 2 ^)A(^iy) 

(- 1 ) 2 e n^n^ n^^^,^^) 
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fe=l fc=l 



using Lemma 3.1 in the next to last equality, has no singularity in the positive 
quadrant {yi > and Zj > with i,j = 1, ...,n}. Indeed we compute the 
residue of this function, near y io = Zj . In the first sum ^ 2E ^ )«, the only 
terms for which the denominator blows up when y io = Zj are those for which 
uu = u>2j ] indeed in the denominator of (|55|), the only factor involving both 
yi and Zj is u 2 j Zj — uji jio y io) which vanishes at yi = Zj if and only if 
u>ii = uj 2jo . 



Setting u 2 j = uiu , the first fraction in fl55|) behaves for y io near Zj , as 



follows: 

A(u l y)A(uj 2 z) 



1 UeKioi^uyt - u Uo y io ) Uio<i(uii yio - uuVt) 

UiKjoi^iZj - u 2jo z jo ) \ [ r ; U:>;.,: ; , - U 2j Zj) 
Ylk&oi^kZk -WlioVio) 

x{an expression not involving Zj ory io } 

1 fUe <i0 (^uyi )Yl i0< i( UJ ii z jo -uj u y £ ) 

uiio{zj -yio)\ Ue^ioi^UoZjo -u u yi) ' 



H i<jo (^2iZi - UJ 2jo Z jo ) H jo<j (^2j Zj - UJ 2j Zj 



+ O(y io -z jo ) 



Hk^j ( U 2kZk - U 2jo Z jo 

x {expression not involving Zj ory io } 

' ' D^i-D^io x { same expression not \ + 0( \ ). 

' * ' lint rr\ hnnrr n i . r\v ry . no o r\ r\-t ro I \ / 



u lto (z jo - y i0 ) K { involving Vio or z jo as above 

The second part in (f)5| ) behaves as: 



fe=i fc=i 



,19 
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k^jo k^i 



n e E"i(hk^) i -("ik!'k) i )*i jj 



i a 



H T 

k=l k=l 



+O(z jo - y 



Jo 

"U JU 



Combining the two estimates, the sum of the terms in Q55D, corresponding 
to u>u = ^jo, can be estimated as follows: 

n n 

_ II II & ^, 3 - A (^2 2)A(a?i y, w 2 z)t 

wi,w a e(C P ) n i=i j=i 
withwi,i =aj 2j0 

wi:=(wii,... ) <2>ii ,...,<<;i 1 n)e(Cp) n 1 = w 2j e Cp 

W 2 : = (w21v.^2j v.^2,n)G(Cp)' I ~ 1 

/ 1 a Wl . o 6 Wl . o J expression independent \ tt tt , nn A 
U-^ ^ X I of y4o and z jo ) 11 °- 11 ^ + 



1 a w 6 w -r-r ^ -r-r f expression independent i 

i,t 11 Wu 11 Wlj I of 77„- and z„- I 



- Vin v ^ 11 11)14 11 u \ of y io and z jo 

= 0(1). 

using condition (|54"D. 

To prove the last statement that lim^^^ E(y, z) exists, note that, in view 
of (4), 

X(y,z) = — + 0(z-y). 
z-y 

Therefore 



aJ)uX{uy,uz) + X (toy , u' 2 
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I h 0(1) for z near y 

/ — ' lo I z — y 

k ^(p ) 

0(1), 



ending the proof of Theorem 6.2. 



7 Fredholm determinants and proof of The- 
orem 0.1 

The next point is to compute Fredholm determinants associated with the 
kernels constructed in Theorems 6.1 and 6.2. First we need a formula gener- 
alizing a classic soliton formula: 

Lemma 7.1. 

1 n 1 n ( \ 



ordered 

Proof: 



<i,j<n 



T 

1 

ordered 

n n 



1 - 

- IT (1 + OiX(y h Zi))T 
T J. J. 

ordered 

-(1+ ^ aiX(yi,Zi)T+ a i a j X(y i ,Zi)X(y j ,z j )T H h ]^[ a» z*)t) 

l<i<n l<«<j<" 1 1 

1+ J2 a i X{VuZt)T + det ' 



X(yj,Zi)T X{ yj ,Zj) 
l<i<n l<«<i<n \ 1 T J T 



H h det(aj )i<ij<n 



, ,/c . X(yi,Zi)r{t) 
aei{0ij + a,j 7— )i<i,j<n- 



T(t) 

This formula is well known in the context of soliton solutions to KdV and 



KP, where it is applied to the r-function tq — 1; see for instance [TO, O, ITS 
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For KdV, we find the customary KdV n-soliton solution, i.e., 

, = 2-^.ogr tt ), 

with t — t + (x, 0, 0, . . .) as in Sect. |] and 

\ 



l<i ,j<n 



T ( t ) = -\ ff e *iX(yu-Vi) r = det^ 7 - ^_ e -E fe: odd* fc (^+^) > | 

to V I V Vi + Vj J 

\ ordered / 

The following is a continuous version of Lemma 7.1, where one replaces 
the determinant appearing in Lemma 7.1 by a Fredholm determinant. This 
statement is due to Adler-Shiota-van Moerbeke |J. 

Theorem 7.2. Given the kernel 

k(y,z) := J D- 1 ((^a^*(t,^)(^ 6 w ,*(t ,«/*))) = ^Y(y,z)r, 

and a subset E = Ui=i[ a 2i-i? a 2i] C R+, the Fredholm determinant of 
k E (y, z) := k(y, z)Ie(z) can be expressed in terms of the original r-function: 

det(J - Xk E ) = I e - A / E d2y (^)r, (56) 
r 

where 

Y(z,z) := limY(y, z). 
If the t- function satisfies for some k > — 1 the relation 

W$t = C kp T, 
then rdet(J — \k E ) satisfies for that same k > —1, 

Proof : We now apply the Fredholm determinant formula 

oo „ „ 

det(/ - \A) = 1 + ^(-A) m ... det (A(z l ,z j )) 1 ^ m dz 1 ...dz m (58) 

m=i ^ - y 

Zl< — <Z m 
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to the kernel A = k E (y,z) := k(y, z)Ie(z), for a given subset E C IR+. It 
suffices to prove for E = [a,b\. At first, the determinant appearing in ( j58|) 
can be computed, according to (|^) and ( |5"T| ) of Theorem 6.1: 



F(z h ...,z m ) := det{k(zi, Zj))i< id < m = - j JJ z*) J r. 

\ ordered / 

We replace the integrals above by Riemann sums, with 

6 — a 

Zi = a + (i — l)fc, = — — -, = a, ^at = A;, 1 < i < N, 



(59) 



so 



oo „ 

det(/-AA; E ) = ^(-A) m / F(z l5 ...,z m ) ^...d^ 

m -Q J a<z 1 <---<z m <b 

N 



m:0 l<i 1 <i 2 <---<i m <N 

We now compute the Riemann sums 

N 

E(- A ) m E (^r 

m:0 l<ji<i 2 <--<im<A r 

2^ (-Adz) , using (53), 

m:0 l<ii<i 2 <---<im<A r ' 



1 N 

— II (1 — A 5z Y(zj, Zj))r, using Lemma 7.1 



ordered 



iTTg-Afay^)^ uging (33) 

Y^ A<5 * y ^V, (60) 



r 

using the fact that 



[y(Zjj 2^'), Y(z m , Zm)] 0, 2j <C 2^,, 
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this is so because 



[X(t,UJZj,UJ 'zj),X(t,u' 'z k ,u;" ' z k )] = 0, 
since 

\ooZj\ <\u) Zf.\ and \u) Zj\ < \u z k \. 
Finally, taking the limit of the exponent in (pO), leads to 

X i 



N ,b 

lim y X5zY(zi, Zi) = X / dzY(z,z) 



establishing (|56|) . This fact together with Theorem 4.2 ends the proof of 
Theorem 7.2. 

Corollary 7.2.1. Given the kernel 

k(z,z') = D' 1 ^?*^, uz)^(t, uj'z') 
for two distinct roots u, uj' of uj p = 1, we have 

det(J - Xk E ) = -e-^fE^x^^'^r. 

T 

Proof : Since uj ^ uj', we have = b w / = 1 and a w > = b w = and thus 
Y) 2^ = 0. Therefore 

z) = X{t, uy, uj'z) and thus Y(z, z) = X(t, ujz, uj'z), 

yielding the result. 

We now define a new kernel 

K(\X'):= l --^% (61) 

where A = 2 P and A' = For any E C R+, let if £ (A, A') = if (A, A')/^'). 
Since 

if (A, AOv^V^A 7 = fc(«, z')\ZdzVdz', 

we have 

det(/-if E ) = det(/-A; sl/P ), 
where i? 1 ^ := {s G R+ | x p G -E 1 }. Thus from Theorem 7.2, we have: 
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Corollary 7.2.2. For the kernel K(X, A') and a subset E C we have 

det(J - u.K E ) = l - {e-vhvv d*Y(z,z)^ T _ (62) 
Moreover, if E is a finite union of intervals: 

r 

E = (J[A 2i _ 1 ,A 2i ], (63) 
i=i 

and if r(t) satisfies for k > — 1, 



w kp r = c kP T > 



E 4 fc+1 ^ + hw£> - c kp ) ) (r det(J - »K E )) = 0. (64) 



then for that same k 

2r 

- A k+1 

dAi 2p 



i=l 



Proof of Theorem 0.1 : The first part of the statement follows from Theorem 
6.2, whereas part (ii) follows from Theorem 7.2. 

8 Virasoro and KP: a new system of PDE's 

To each k G Z, we associate the weight k and to a constant the weight 
0; thus 

/ m \ m 

weight [cJljV ) =$>• 

Given the operator 



i=i 



Q - E a kl---krn J kl J S ■■■ J kl 

we define its "principal symbol" 

cr(Q) = sum of terms in Q of maximal weight 

and the "highest gap" 7q is the difference between the maximal weight and 
the next to maximal weight in Q. 
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In the next lemma, we consider t-operators of weight at most k of the 
form 

Qk = ^akiJ^ + ^2 h,i+f + d k wi th a kk = 1, (65) 

i<k i+j<k 

with a constant d k , and differential operators A k in parameters a of precise 
weight k. We also need to give a precise meaning to 

1 1 

Pi(d A ) := Pi(A u -A 2 , -A 3 , ...); (66) 

the commutative variables ^ get replaced by the non-commutative operators 
\Ai with the following ordering: if i ± < i 2 < ... < it, then 

1 A 

ti 1 ti 2 ...ti e I > - —J\i t ...J\i x . 

Lemma 8.1. Let a KdV r-function r(t,a), depending on t e C°° and an 
extra-set of parameters a, satisfy equations of the type 

Q kT = A k r, with k = 1, 2, (67) 

with Q k as in (65), then the following holds : 

ifb k)k = for all k = 1,2,..., then r(0,a) satisfies a hierarchy of bilinear 
PDE's, of the form 

T ■ A n AlT - A n T ■ AlT - ^ Pi0a)t ■ Pj(d A )T 

i+j=n+l 

+ (terms of weight < n + 1 — 7^ = (68) 

for n = 3, 4, 5, .... 

If one shifts t\ \— > t± + x in the equations (67), then x figures in the lower- 
weight terms of (68). 

if some b k>k ^ for some k = 1,2, and, provided some determinant in- 
volving the b kjk is 7^ 0, then r(0, a) satisfies a similar hierarchy of PDE's: 

r ■ B nl T - B n r ■ B\T — ^ Pi0B)r ■ Pj(d B )r 

i+j=n+l 

+ (terms of weight < n + 1 — 7^ = (69) 

for n = 3, 4, 5, 
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where p n {ds), £>« arid B n -i t i are appropriate linear combinations of 
A je ...A jl , for < ji < j 2 < j' a = n, 



and 



7 = mm7 Q 

kez Vfc 



Proof: Consider the infinite set 



<Sk = {^i > k 2 > ... > kt with ki E Z and k a = k for arbitrary £} 



and the finite set P fc of partitions, defined by 

t 

S k D Vk = {k\ > k 2 > ... > 0, with hi E Z and k a = k for arbitrary £} 

l 

= {Young diagrams Y of weight fc}. 
In terms of the "principal symbol" a(Q k ) of Q k : 

a(Q k ) = jl 1) + b kk Y^-.jWj™:, k>l, 

i+j=k 

we have for given Y = (hi > k 2 > ... > 0) e V k 

Qkf-Qk! = <?{Qk l )~-<?{Qk 1 ) + lower weight terms (70) 

(il>...>im)&S k 

with the c^'"^ being functions of the fefc^ only. For a given Young diagram 
Y = (ki > k 2 > ...), we define 



t(i) ._ t(1) tW. 
J Y ■— J k e ■■■ J k 1 i 



with this notation 



Qki---Qki 



Jy^ + E c y ' + l° wer weight terms + c. (71) 



t=o 

Y'ev k 



35 



Adler-Van Moerbeke : Random matrix February 9, 2008 #1 §8, p. 36 



Thus, for given k > 0, we get a system of (^Pfc)-equations in the unknowns 
Jy \ y G Vk, which is invertible, provided 

det (l+(4') Y >,YeP k ) t^O. 
Hence, for each Y = (ii > i 2 > ...) G Pfc the equations (71) yield 



<>1 ("m 



(h-h)eV k 



t=o 



K 1 ) 



+ lower weight terms in 
with fc > 

We now proceed by induction: for k — 1 i.e., V\ = {(1)}, we have 
Qi = + (lower weight terms) 

and thus 



d 



+ constant. 



t=o 



Therefore for each Y = (i\ > i 2 > ... > i m ) G Vk, k > 1, 



j 1 >...>j e >i 
i<j 1 +...+j e <k 



+ constant. 



t=o 



Therefore for a r-function satisfying (67), we have 



t=o 



ii>. ..>«>i 
l<J'l+...+Jf<fc 



t=0 



+ constant, (72) 



since, for instance, 

Qj2Qji T = Qj2-Aji T = J ^ i jiQj2 T = •Aj 1 A.j 2 T. 
Note r also satisfies the .fTP hierarchy 



i+j=n+l 



(73) 
n = 3,4,5,..., 
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whose terms all have the same weight n + 1. 

For instance, the expression p n (d) in the KP-equation, namely 



Pr, 



d_ ld_ 



(jL 1) ) ka ...(j[ 1) ) ki 



^ /•',.:.../•' j !f I '' ' 

ki > 



Ifei' 



leads to a contribution 



£;„!.. .LI — 



E 



E 



j 1 >...>j e >l i Ki ^v 
Y.ih= n \T." ik i= r 



1-31 

a,...,a,a—l,...,a—l, 

Hi r kr k a \...ki\ 



+ lower weight terms 
\ 



A.j 1 ...A._ 



■jt 



J 



+ lower weight terms. 



When bk t k = 0, for all k = 1, 2, then the proof is much simpler; indeed, 
then formula (70) becomes 



Qkv-Qki = ■■■Jkx + l° wer weight terms in jjp, 



(74) 



which leads at once to 



J (1) J (1) r 



= Akf-AkiT +( lower weight terms in Ak) 
t=o t=o 



Upon shifting ti*—*ti + x, the variable x will only figure in the lower-weight 
terms of (65), using = 0, and thus x appears in the lower-weight terms 
of (68) only. This ends the proof of Lemma 8.1. 



In the proof of Lemma 8.1, one needs to compute expressions like (|70|): 
for (ki > ... > k e ) G V k 



Qki • --Qki 



t=0 



J (1) J (1) + 

J k e ■■■ J k 1 + 



E 

(ii>->im)£T>k 
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i.e., typically, by virtue of the form (65) of Qt, one is dealing with products 
of normal ordered sums 



(}2 * r : ) . . .( ^ : 

\i+j=k e / \i+j=ki 



which are then evaluated at t — 0. They can actually be computed in terms 
of graphs, as explained in Lemma 8.2. 

Consider the following generating function 

T(z) := : jW(z)jW(z): , with J^(z) = E J^z"*" 1 and = /i 



n£Z \i+j=n / 

= E^~ n ~ 2 - 



Upon shifting t n i— > t n + /x n /n 



j«w=x;^- i - i +E^" 1 



and thus 



t=o 



iez 

i>l 



Also 



t=Q 



Consider the set of graphs 



where is a graph with n numbered vertices subjected to the rule: out 
of every vertex, you have at most 2 links leading to another vertex. Each 
graph = lj i G^ consists of a number of connected components as in 
the figure below 



p(n) = 
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We will be using the following series: 

1 



E 



iz 1 1 w t 1 = 



i>0 

Lemma 8.2. 



{z — w)' 



and 



Of - 1) _j_ 2 j_ 2 



i>0 



— u>) 



4 ' 



no-no = E IU 

Gi n) erw i 

where each G^} contributes a factor in each term: 



G (n): 



if = {k} : 

n (n) _ tl ll 

U ai - ' ' ■ 

n (n) _ m-i m 2 

^ai - ' = ' : 



fa ai = T{z k ) 

let p rCll/ \ t("1W \ r-ra-1 1 



/ 0a4 = jW(^ 1 )J«(^)IKi 1 



K-2£ i+1 ) 2 



/Oai 2( Zm ,-2 m2 ) 4 



Proof : The proof proceeds by induction on m, the relation being obviously 
true for m — 1. 



Corollary 8.2.1. We .have 

T(zi)...T(^ 



t=o 



G w er („) 



t=o 



where 



if G 



(n) 
(n) 



{*} = 



<3 



mi m-z 



r>\ n > — 

Mn) _ 



9 Proof of Theorem 0.2 

Remember from section that 



tot 



f r (n) 


t=0 


f r (") 


t=0 


oti 


t=0 



■^(^fc) 2 



^(-^j^co nr=i 1 ( Z£ .-i. +l )2 



t=0 2(z mi 2m 2 ) 4 



^ = i, ""'-^ ■ v\*)) ■ E r * z 



2i 



"dz 



i>l 
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with 

Viz) = - z +-z 3 ^0,m = deg V = or 2. 
2 6 

Lemma 9.1. Consider a function ^f(x, z) satisfying conditions (12) and (13) 
of Theorem 0.2; i.e., ^(x,z) is a solution of the linear partial differential 
equation 

A g V{x,z)=A x V(x,z) (75) 
with holomorphic \I/(0, z) in z 1 , subjected to 

(aA 2 + bA + c)^(0,z) = z 2 *(0,z), with tf(0, z) = 1 + 0(^ 1 ). (76) 

Then ty(x, z) specifies an infinite- dimensional plane 

W:= span c {(J-J *(x,z)\ x=0 , j = 0, 1, 2, . . .} G Gr°(x), (77) 

invariant under z 2 and A, 

z 2 W C W and AW C W, 

such that 

W x := e~ xz W = span {e~ xz (J^ z),j = 0, 1, 2, . . . for hxed x). 

with *$?(x, z) = e xz {l + 0(2 -1 )). Also, there exists a potential q(x) such that 

Remark : Another way of stating Lemma 9.1 goes as follows : the tau- 
function r(t) associated with 

W x,t = e -xz-YZti* w ^ with w ag in (j 7 ^ 
and the corresponding wave function 



V(x,t,z) = e xz+ ^ u 



;r{t-[z-'\) 



r(t) 

are such that ^f(x, 0, z) = "$?(x, z). 

i5^< r (o) j g ^.^g ce ^j Q £ ^Yi e Grassmannian Gr. 
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Proof : First observe that (m = or 2): 

Id a 6 

A z = -z~ m+l — + c m z + 0(1), with c = - and c 2 = ^ (78) 
A Oz 4 4 

with 

[A z ,z 2 }=z 2 -™. (79) 

Consider the function \I/(x, z), solution of (75) and (76) ; then the plane 
W e Gr, generated by *(0, z), 

W := span {A k 2 V(0,z), k = 0,1,2,...}, (80) 

belongs to the big cell Gr^°\ because 

A k ^(0,z) = (c m z) k (l + O(z~ 1 )). 

The plan W has the following property 

A Z W C W and z 2 W C W. 

The first inclusion follows at one from the construction (80) of the plane W, 
while the second requires some argument. It suffices to prove 

fc+2 

z 2 A k ^{0, z) = J2 CiA l V{0, z), (81) 

i=0 

which is done by induction. The case k = is merely condition (76) and we 
assume (81) up to including k; we compute 

z 2 A k+1 ^(0,z) 

= A k+1 z 2 ^(0, z) + [z 2 , A k+1 ]^(0, z) 

= A k+1 z 2 ^(0, z) + Y!U A'[z 2 , A}A k ~^(0, z) 

= (aA k+3 + bA k+2 + cA k+1 )^(0, z) - £- = o A's 2 -"^*-** (0, z) 

using (76) and (79). 

= YaIq C i Al ^(^^ Z )^ for SOme C i- 
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The process of letting W evolve according to the i^P-vector fields defines a 
T-function and a wave function ^(x, t,z) = e xz+ ^ tiZ% (1 + 0(z~ r )), in terms 
of which the same plane W can be expressed as 

W= span{Q^ V(x,t,z), j = 0, 1, 2, . . . for fixed (x, *)}. (82) 
Therefore, whenever AVF C W for a differential operator A = A z , we have 

Att(a;,f,z) = X^Gr) *0M,4 (83) 



Also, remember, to any operator A G C[z,z we associate a pseudo- 

differential operator P4 



— oo<i<oo — oo<i<oo 

0<j<OO 0<J<OO 

such that 

A*(a;, t, 2) = P A ^(x, t, z). (84) 

Upon comparing (83) and (84), we have that Pa is a differential operator for 
that i.e., P4 = (Pa)+. 

Therefore, in view of the explicit form (9) of A and the representation 
(23) of M, we have 

P A = \{M + V'(L))L- m+l + E,>i c. 2l L~ 21 

= {\ML- m ^ + c m L) + (85) 
= (S{\xD-™^ + c m D)S- 1 ) + + ± E fe >i kt k (L k ~ m ) + 

and 

P Z 2 = L 2 = (L 2 ) + = D 2 + q(x,t). (86) 
Setting t = in (85) leads to 



A\t=Q 



= (s(^xD- m+1 + c m D)S-^ 
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= ((1 + a^D- 1 + -..)(\ xD- m+1 + c m D)(l - a^D- 1 + . . .))+ 

= jS mt2 D + ^(x+^)6 m , D = A x . (87) 

Therefore, combining (84) and (85), fy(x,0,z) is a solution of the first order 
partial differential equation 

A g V(x,0,z) = A x V{x,0,z) (88) 

with 

tt(0,0,z)=tt(0,z), 

where ^(0,z) satisfies (76). But the function ty(x,z), alluded to in the 
statement of Lemma 9.1, satisfies, by assumption, the same equation (75), 
with the same initial condition (76). The solution of such a linear partial 
differential equation is unique, because the relevant Jacobian does not vanish, 
namely the coefficient of D in A x . Therefore we conclude ^(x, 0, z) = ^(x, z) 
and the plane (77) generated by ^(x, z) has the properties stated in Lemma 
9.1. 

Lemma 9.2. A function *f?(x,z) satisfying conditions (12) and (13) of The- 
orem 0.2, defines a unique KdVr-function r(t) = t(x + t 1 ,t 2 , . . .) satisfying 
the Virasoro constraints 

^(W^ + aW^ + (3Wi%) = c k r for even k > -m, (89) 

for some constants c k ■ 

Proof : From the proof of Lemma 9.1, we have that 

P Z 2 = L 2 and P A = \{ML~ m+1 + \aL~ m+1 + ^/3L _ro+3 ) 

Zi Zi Zi 

are differential operators and thus (P A L k+m )_ = for all even k > —m. 
Then from the ASV-formula (see (48) in lemma 5.2): 

— (M n L n+£ )_^ ,__ n ^W? +1 t 



= (e-"-l)- 
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it follows that for all even k > —m 

Q (P A L k+m )-V 

1 (ML k+1 + f L k+1 + f L fc+3 )_^ 

and thus (89) holds, ending the proof of Lemma 9.2. 

Proof of Theorem 0.2 : Consider now the shift 

a (3 
si = ti + -, 3s 3 = 3t 3 + -, all other Sj = t i? 

and define 

f(s) = r(f). 

For even /c > — m, we have 

„2 



— 2 

where is H4, with t replaced by s. Note the term ^-5k-2 comes from 



W ( V= S j + ... = (t 1 + a/2) 2 + ... = w[ 2 l + aW { _}l + ^ + .... 
Therefore the equations (89) in Lemma 9.2 get converted into 

W^f = c k f, for even k > —m, (91) 

with different constants c^. 

Let i&(x,s,z) and ^*(x,s,z) be the wave functions defined by the r- 
function f(s); note for KdV, we have ^*(x, s,z) = ^f(x,s,—z). Using this 
fact, combined with the requirement a + b + = a_5_, implies the proportional- 
ity of the functions 

^2 a^*(x,s,ujy) = ^ a u V(x,s,-wy), 
^e{+,-} u>e{+-} 

b^(x,s,u'y) = Yl b^(x,s,uj'y)=:<S>(x,s,y). 

^'G{+,-} «'€{+,-} 
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The kernel 

flows off the kernel (15) at s = (|, 0, |, 0, 0, . . .). According to Theorem 0.1, 
using the Virasoro equations (91) for f , the expression 

f E ■= fdet(7 - \K E S ) = e - x ^ dz Y ^ z ^f(s) (92) 

satisfies the equation 

2r k r, \ 

-2 ]T A? + + -(Wf > - c fc ) f £ = 0, for even fc > -m, 
i=i ' / 

and, upon undoing the shift above, r E : = rdet(J — \K E t ) satisfies 

(- 2 E + 5WP + + wSi - *)) - £ = 0, (93) 

for even fc > —to, 

with principal symbol 

<r(wf > + aWt% + (3W^ 3 - c k ) = (3W^ 3 if f3 ? fc even > -2 

= a^J if /? = fc even > 0. 

We now apply Lemma 8.1 to the case where all b^k = 0. Since we are 
in the KdV situation, we only need for odd kf, this is obtained 

upon considering the constraints QkT E = AkJ E for odd k > 1 and the KdV 
equations (73) for odd n > 3. So Lemma 8.1 is applicable and r E satisfies 
the hierarchy (68) of KP equations for odd n > 3, with the commutative 
t-partials replaced by the non-commutative ^-operators. 

Note these equations do not have constant terms. Indeed, shifting the s- 
variables back to the t-variables in (92) leads to the appearance of exponential 
factors in (92); differentiating that new expression (92) with regard to the 
boundary points Ai and letting all Ai tend to oo implies that in the limit all 
terms in (68) tend to zero. 

Finally, because of the fact that in r E = rdet(J — \K E X ) the factor r(t) 
is independent of the boundary points Ai and that the equations are free of 
constant terms, the Fredholm determinant itself is a solution of equations 
(16). This ends the proof of Theorem 0.2. 
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10 Examples 
Example 1 (Airy kernel) 

Given an entire function U(y) growing sufficiently fast at oo, consider its 
"Fourier" transform, 

/oo 
e -U(y)+yu dy . ( Q4 ) 
-oo 

define the associated function p{z) and the differential operator A := A z : 
p(z) = -^e^^'^v^ and A z := p{ z )-^.^-p{ z )-\ (95) 



Proposition 10.1. The function 

z) := p(z)F(x + U'{z)) (96) 
satisfies the following three equations: 

A z ^(x,z) = ^(x,z), (97) 



(u'(-^) - x^j *(x, z) = U'(z)V(x, z). 



(98) 



and 

U'(A)V(0, z) = U'(z)V(0, z), (99) 
with the following asymptotics 

#(0,z) = 1 + 0(z- 1 ) for z / oo. (100) 

Proof : At first note that 

z) = P {z)^F{x + U'(z)) = P(z)jJ^Y/(x + U'(z)) = B*(x, z) 

Integrating -jL e ~ u ^ +yu over R shows that F(u) satisfies the following differ- 
ential equation 

U'(J-^F(u)=uF(u). (101) 
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Setting u = x + U'(z) in this equation and multiplying with p(z) yield: 



The method of stationary phase applied to the integral (94) and its deriva- 
tives leads to the following estimate, upon Taylor expanding U(x) around 
x — z, and upon using (95): 



thus ending the proof of proposition 10.1. 

Theorem 10.2. ForU(z) = z 3 /3, the function fy(x,z), defined by (96), sat- 
isfies the conditions of Theorem 0.2 and the Fredholm determinant f(Ai, A 2r ) 
:= det(I — AfCf ) satisfies the hierarchy of bilinear partial differential equa- 
tions in the for odd n > 3: 






p(z)-V(l + 0(l/z)), 



(102) 



f ■ A n AJ - A n f ■ AJ - PMf-Pj(-A)f = 0, 



(103) 



i+j=n+l 



with 




n 



— 1, 3, 5, 



(104) 
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the variables appearing in the Schur polynomials pi are non-commutative 
and are written according to a definite order. Finally, the first equation in 
the hierarchy (16) takes on the following form : 

Al - 4(^3 ~^ R + 6(Ai?) 2 = 0, (105) 
for 

2r 



R := Ax log / = ^ — 



When E = (— oo, A), the function R = A\ log/ = ^ logdet(7 — K E ) satisfies 
R'" - AAR' + 2R + QR' 2 = (Painleve II) 

Proof : The Painleve II equation for the logarithmic derivative has been ob- 
tained previously by Tracy and Widom |28]] ; equation (105), which leads to 
Painleve, is new. For U(z) = |z 3 , we have a = 0, (3 = 4, m = 2, and 
V(z) = §z 3 ; also 

A:= — (^- + 2z 2 ) --z- 2 and A= (106) 
2z \dz J 4 dx K ' 

Then in terms of the Airy function 

/OO 3 
e- v - +yu dy, 
-oo 

ty(x, z) has, according to proposition 10.1, the following expression 

^(x, z) = -^-=e~i zi 'y/zF(x + z 2 ) 

'71 



and is a solution of 

A$(x, z) = A$(x, z), 

with ^(0,z) satisfying 

B 2 ^(0, z) = z 2 ^(0, z) and ^(0, z) = 1 + 0{z~ 1 ) 
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Setting b + — 1 and 6_ = 0, we find for u) and K E (y,z) in (15), 



$(x,u) = ^-=A(x + u) 



i r 

K E {y,z) = I E (z)— A(x + y)A(x + z)dx. 
2n J 

Then f(Ai, . . . , Az r ) '■= det(J — \K E ) satisfies the hierarchy of equations 
(16), with lower weight terms; the A n are defined by (107). However, upon 
rewriting the variables of p n in an appropriate way, all lower weight terms 
can be removed, as follows from a combinatorial argument. 

Finally r E := rdet(J — XK E t ) satisfies equation (93) for a = 0, (3 — 4 
and m = 2; in particular for k = 2 and 0, one finds, upon dividing by r: 

A 1 ]ogT» = Z I > >^ + 2-J-ji ogr £ + ! 



^3 logr 1 




1 

16 



from which the partial derivatives log r J , ^ log r J and log r J at t = 
can be extracted. Putting these partials in the KP-equation (73) leads to 
the equation (105). In the particular case of a semi-infinite interval, one finds 
the Painleve II equation, as stated in Theorem 10.1. 

Example 2 (Bessel Kernel) 

In the next theorem, we investigate the PDE's for the Fredholm determinant 
of the Bessel kernel, depending on a parameter u; special values of v yield 
the sine kernel, which was originally investigated by Jimbo, Miwa, Mori and 
Sato [T3|. The equation (109) is new, although in the case of an interval the 



Painleve equations where obtained in [jTjJ and [28[ . 
Theorem 10.3. Given the kernel 



i r 1 

K E (y,z) = --I E (z) sJ v {s^)J v {s^)ds, (107) 
* Jo 
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the Fredholm determinant f(Ai, Ajr) '■= det(J — \K E ) satisfies the hier- 
archy (16) of bilinear partial differential equations in the Ai for odd n > 3, 
with 

^ , n+1 Q 

An^TA 2 ^7-, n= 1,3,5,...; (108) 
i=i OAi 

The first equation in the hierarchy (16) for F := logdet(J — \K E ) takes on 
the following form: 

Ai - 2A\ + (1 - v 2 )A\ + A 3 - F-A{A!F)(AlF)+6{AlF) 2 = 

(109) 

when E = (0, A), we have for R := —AiF = -A ^ logdet(J - XK E ), the 
equation: 

A 2 R"' + AR" + {A - v 2 )R' -- + ARR' - 6AR' 2 = (Painleve V) 

2 

Proof : Pick V(z) = —z; then m = 0, a = —2, f3 — and 

Id -Id 
A z = -z(- 1) and A x = -(x-l) — 

2 K dz J 2 v J dx 

We look for a function ^(x, z) satisfying 

A z V(x,z) = A x V{x,z) (110) 
with initial condition \l/(0,z) satisfying 

(AA 2 z -2A z -u 2 + h^(0,z) =z 2 V(0,z), *(0, z) = 1 + O(-) (111) 

4 ^ 

The solution to the differential equation (111) is given by 

V(0,z) = B(z) = e^H u (iz) 

g 22^+l/2 roo z -v+l/2 e -uz 



v(-v + 1/2) j 1 {u 2 - iy +l / 2 

- 1 + 0(|) 
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with e = i^/Vfee™^ ,-\<v<\. Then 

V(x,z) = e xz B((l-x)z) 

satisfies (110); from Theorem 0.2, ty(x,z) satisfies a second order spectral 
problem, which can explicitly be computed : 



d 2 (v z - 



2 — i x 

4' 



^(x,z) = z 2 ^(x,z). 



^dx 2 (x — ly 

Picking b + = e~ im/ ^ /2y/n and 6_ = ib+, yield for (15): 

*(x,z) = —^e- z q(x,z) + —^e z q(x,-z) 
[ ^^U(l-x)iz), 



and 



K, = I E (z) I 2 yV4^i/4 dx 



x 



I r(l-x)i 

= --I E (z) / sJ v (sy/y)J v (s\fz)ds. 
* Jo 

The special value x — i + 1 leads to the standard Bessel kernel: 

i r 1 

Kf +l = --Ie(z) J sJ u {s^)J v {sVz)ds 

~ Ie{z) W^v) 

From Theorem 0.2, the Fredholm determinant 

f(A±, A 2r ) :=det(/-AXf +i ) 
satisfies equation (16), with 

2r 

n+1 

A. 2 



*— nil d 

• A « = 22 A i 2 o77> n = 1,3,5, 



i=i 



51 



Adler-Van Moerbeke : Random matrix February 9, 2008 #1 §10, p. 52 



Picking the special value x — i e 1, leads to the shift t\ \— > t\ + i + 1; setting 
a = — 2 and /3 = in (93) and shifting ti, as above, leads to 



\i>l 



a ia 2 , 9 \ , ,, 1 / d 



Alogr- = il > + + v'-l^ J logr- + j logr* ; , 

expressing the partial derivatives appearing in the first KP-equation for p = 2 
at t = in terms of the operators „4i = YlT A^- and .4.3 = YlT M" leads 
to the partial differential equation (109), which for E = (0,A) leads to the 
Painleve V equation, ending the proof of Theorem 10.3. 

Sine kernel : Note that for v = ±1/2, 

2 k ~(V, Z ) ^±1/2/ 2 2n 



7T ^/p 



1 



x sin sin £Z 
cos COS , 

eta; 



7T Jo y 1 / 2 z 1 l 2 

— -— / ( e -«» Te ^)(e« Te - ,K )(li 

1 . sin x(y — z) sin x(y + 2) . 
2n y — z ' y + z 

Note this kernel can be obtained in a direct way, as follows: assume r = 1, 
p = 2, ai = 61 = 1 and a_i = 6_i = ±1 in (5); one finds 

kt(y,z) := fc t (y,z) 

= / (V*(ix,t,y)±V*(ix,t,-y))(ty(ix,t,z)±V(ix,t,-z))dx 
Jo 

where 

= e ixz e^ tkzk and = e ~ lxz tkzk . 
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Note the kernel k t (y, z) flows off the kernel at t — 0: 

PX 

k± Q (y, z) = {e~ lxy ± e lxy ) {e lxz ± e~ lxz ) dx 
Jo 

^ / sin x(y — z) ^ sin x(y + z) \ 
V y-z " y + z ) ' 

Moreover 

det (/-AAi7 al )=det (/ - 2\k^f) . 

It satisfies 

|^_ a 2fc+i + I jP)^ det ^ _ AA;l- a ' a] ) =0 for fc > 0, (112) 

since J^l = for k > 0. Also for K(X, A') = (l/2)fc(^, ^O/v 7 ^, A = 
A' = z' 2 , we have 

(-^ +1 ^ + ^2 ) )det(/-2A4r)=0. 
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